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Abstract 

C*") , In this paper we review two concepts directly related to the Lax representations: Darboux trans- 

formations and Recursion operators for integrable systems. We then present an extensive list of 
integrable differential-difference equations together with their Hamiltonian structures, recursion oper- 
ators, nontrivial generalised symmetries and Darboux-Lax representations. The new results include 
multi-Hamiltonian structures and recursion operators for integrable Volterra type equations, integrable 
discretization of derivative nonlinear Schrodinger equations such as the Kaup-Newell lattice, the Chen- 
Lee-Liu lattice and the Ablowitz-Ramani-Segur (Gerdjikov-Ivanov) lattice. We also compute the weakly 
nonlocal inverse recursion operators. 

<n : 

1 Introduction 

GO ' 

The aim of this paper is to give a comprehensive account of multi-Hamiltonian structures, recursion op- 
erators and Darboux-Lax representations for a wide class of integrable differential-difference equations. 
Some of these results are well known, but scattered in literature. In many cases we have completed the 
picture providing explicit expressions for Hamiltonian, symplectic, recursion operators and Darboux- 
Lax representations. The Lax representations of nonlinear differential and difference equations play a 
central role in the theory of integrable systems. It enables one to apply the inverse transform method for 
construction of exact solutions and study the asymptotic of the initial value problem. It also enables one 
, to construct a recursion operator which generates a infinite hierarchies of symmetries and conservation 
laws. Currently for a given equation there is not any general method to find its Lax representation. The 
most successful approach is the Wahlquist-Estabrook prolongation procedure [1J. Recently, Mikhailov 
' and his coauthors have tackled this problem from a different angle |2H3l[4j[5). They studied possible 
, reductions of a general Lax representation using the reduction group approach (6) |7J [8) and further 
leading to a classification of the Lax representations and corresponding integrable equations. 
The concept of Darboux transformations is originated from the classical differential geometry [9J. Appli- 
cation of Darboux transformations to the corresponding Lax representation leads to Backhand transfor- 
, mations and to a generation of new exact solutions for the integrable system. Backlund transformations 
can be regarded as an integrable system of differential-difference equations on their own rights. These 
differential-difference equations play role of infinitesimal symmetries for the integrable partial differ- 
ence equations which can be obtained from the condition of the Bianchi commutativity of the Darboux 
transformations. 

Evolutionary differential-difference systems are the main object of our study. Our notations are standard. 
We illustrate them on the well-known example of the Volterra equation HUI (cf. Section [47T1 for more 
algebraic properties of this equation), which we can write in the form 

ll t = U(ll\ - U-\). (1) 

Here the dependent variable u is assumed to be a function u(n,t) of a lattice variable n e Z and a 
continuous variable f and we denote 

n t = d t u{n, t), u j = u(n + j, t). 
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We shall omit the subscript index zero and use u instead of uq. The Volterra equation |Q} encodes an 
infinite sequence of differential equations 

d t u(n, t) = u(n, t)(u(n + l,t)- u(n - 1, f)), «eZ. 

Equation (TJJ possesses an infinite hierarchy of symmetries, or in other words, it is compatible with an 
infinite sequence of evolutionary equations of the form 

U tm = KmiUm, Um-i, U\- m , U- m ), m € N 

where t\ = t, K\= u{ii\ - u_i), 

K2 = U (u\U2 + ll\ + UU\ - MM_i - M^j - U-\ll-2^ 

and K m are certain polynomials of variables u m , U m -\, . . . , U\- m , U- m . The compatibility means d tm {d t u) = 
d t (d tm u) or the vanishing of the Lie bracket [K\, K m ], which is defined as 

[K\,K m ] := K mi ,(K\) - Kn,{K m ) = 

where K m + denotes the Frechet derivative 




and S denotes the shift operator, such that S(u k ) = u k+ j and for any function f(u kl , . . . , u kl ) we have 
S(f(u ku . . . , UjtJ) = f{u h+j , u h+j ). 

Symmetries of the Volterra equation can be generated by the recursion operator "R 

n = uS + u + u x + uS' 1 + m(ui - U-x){S - i, (2) 

namely 

K m +i = ^'"(Ki) 

The vector fields corresponding to the Volterra equation and its symmetries are difference polynomials, 
i.e. elements of the difference polynomial ring R = [C, u,S] which is a multi-variate ring over the 
field of complex numbers C with infinite number of indeterminates u k , k e Z and equipped with the 
automorphism S. The corresponding field of fractions T — (C, u, S) is a difference field of rational 
functions of indeterminates u k , k e Z over C, which inherits the automorphism S. Difference operators 
are defined as finite sums of the form £L a k S k , where a k e f (the above defined Frechet derivative K m + is 
an example of a difference operator). Difference operators act naturally on elements of T ■ Similar to the 
differential case the elements of the ring R, field F and difference operators are called local polynomials, 
local functions and local operators respectively. 

The recursion operator fi given by 10 is not a local operator. It contains the local part and the term 
u{u\ — U-\)(S - l)~ l u~ l . Action of H can be defined on those elements of !7~ which belong to the the image 
of the difference operator u(S - 1). We can directly check that K\, K2 e Im (u(S - 1)). Similar to IfTTH it 
can be shown that all K m+ \ - ^'"(Ki), m e IN are difference polynomials. The action of K is not uniquely 
defined on the elements of Im (u(S — 1)). Indeed, the base field C is the kernel space of the difference 
operator n(S - 1) and thus for a = u(S - l)(b), b e T we have (S - l) -1 (M -1 a) = b + a, where a e C is an 
arbitrary constant. Obviously, acting by the recursion operator K defined by 10 on K\ we get 

7?(m(«1 - M-l)) = U (ll\U2 + u\ + UU\ - UU-i - ll 2 _ x - U-1II-2} + au(ll\ - U-\). 

Action of the recursion operator H is well defined on the sequence of the quotient linear spaces 

TCi = Span c (Ki), % n = SpanJKi, K 2 , K„,)/Span c (K 1 , K 2 , • ■ ■ , K m -i) 

and in what follows describing the result of the action of a recursion operator on a symmetry we shall 
give one representative from the corresponding co-set. 
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The recursion operator K given by 10 is a pseudo-difference operator. It can be represented in the form 
% = BA~ l , where A and B are difference operators. For example one can take A = *H\, B = 'H2, where 
"Hi and ^2 are two Hamiltonian operators for the Volterra equation (see Section 4.1). Notice that the 
pseudo-difference operator H is a sum of a local (difference) operator and a nonlocal term of the form 
P(S — 1) _1 Q. We say that a pseudo-difference operator is weakly nonlocal if it can be represented in the 
form 



where A, P, and Q, are difference operators and the sum is finite (a similar terminology has been first 
introduced in the study of pseudo-differential Hamiltonian operators llTZl ). Thus the recursion operator 
H given by ((J) is a weakly nonlocal pseudo-difference operator. Moreover, it is easy to show that 
K m , m e N is a pseudo-difference operator. In the majority of cases studied in the present paper the 
recursion operators are weakly nonlocal. Exceptions include the recursion operator for the Narita-Itoh- 
Bogoyavlensky lattice [13] (see Section 4.5). 

We notice that surprisingly, in the case of multi-component systems of integrable difference equations 
with weakly nonlocal recursion operators often the inverse recursion operator is also weakly nonlocal. It 
enables us to generate infinitely many local symmetries corresponding to the inverse flows. For example 
the Heisenberg ferromagnet lattice Ifl4ll (cf. Section l4.17l for more algebraic properties of this equation) 

Ut = (U - V)(U - Ml)(«l - v)~ x 
Cf = (U - v)(V-i - V)(U - V-i)~ l 



possesses a recursion operator 

( (u-vf 
<fc — (ui-v) 2 



»_ 2(u-h 1 )(p-p_ 1 ) 

(tt-0_l)(ttl-») 

(v-v-lf 

{u-v-if 



+ 2K m {S-l)- 1 Q 



(1) 



where 



(»— V)(U— Mi) N 
U\— V 

(u-v)(v-i-v) 



V—V-l 



11— III 



v)(u-v-i)' (u-v)(ui-v) 



The operator *R is weakly nonlocal and it has a weakly nonlocal inverse 



I (u-vf 



( U - Vl y 



(il-»-l) 2 

(11-1 -vf 

(U-V) 2 Q _ 2(»-lLi)(g-Di) 
(ll-V-i) 1 (u-j>i)(u_i-o) J 



2X ( - 1) (vS-l)" 1 Q l 



(-1) 



where 



f (u-v)(u-i -u) y 

H_i —V 

(u-v)(v-vi) 

U-Vi 



Q 



(-1) 



( 



V—Vy 



(U-V)(U-V{)' (; 



U-U-i \ 

n-v)(u_i-v) J ' 



Thus, the Heisenberg ferromagnet lattice has infinitely many local symmetries ^'(X' 1 ') and ^~'(X^ _1 ') for 
all / e INT. The phenomenon has been explored for the Ablowitz-Ladik lattice and the Bruschi-Ragnisco 
lattice in fl5l . In this paper, we compute the weakly nonlocal inverse recursion operators for all multi- 
component integrable differential-difference equations if existing. Such inverses do not exist for scalar 
nonlinear integrable differential-difference equations. For a given weakly nonlocal difference operator, 
to answer whether there exists a weakly nonlocal inverse operator is still an open problem. 
The arrangement of this paper is as follows: First we review two closely related topics concerning the 
Lax representations. One is the Darboux transformations of the Lax representation, from which we 
derive the integrable differential-difference equations. Another topic is to derive the recursion operator 
for the resulting equations using the Darboux transformation. We illustrate the methods through two 
typical examples: the well-known nonlinear Schrodinger equation and a deformation of the derivative 
nonlinear Schrodinger equation corresponding to the dihedral reduction group D2. 
We complete the paper with a long list of integrable differential-difference equations, where we list 
equations themselves, their Hamiltonian structures, recursion operators, nontrivial generalised symme- 
tries and their Lax representations. We also include partial results on their master symmetries. For 
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some equations we add further notes concerning the links with other known equations and the weakly 
nonlocal inverses of recursion operators if existing. The list is far from being complete. 
We mainly refer to the source where we learned about each system although some attempt has been 
made to track the original contributions. Compiling the list we have verified and made consistent the 
objects collected from the vast literature. Our list also includes a number of new results (to the best of 
our knowledge): 

• The Hamiltonian operators, symplectic operators and recursion operators given in section l4~4l for 
equations ll68t-ll70t, and the relations among them; 

• The Hamiltonian operators, symplectic operators and recursion operators for the Kaup-Newell 
lattice (section 14.14b . the Chen-Lee-Liu lattice (section 14.151 and the Ablowitz-Ramani-Segur 
(Gerdjikov-Ivanov) lattice (section [4.16b ; 

• All weakly nonlocal inverse recursion operators if existing (except for the cases of Ablowitz-Ladik 
lattice and Bruschi-Ragnisco lattice, which were known). 

2 Lax representations and Darboux matrices 

With an evolutionary nonlinear partial differential system 

u t = F(u,u x ,...u t .... v ), ueC m (3) 

solvable via the spectral transform method (1611171 [T51 one associates a pair of linear operators 

L = D x - U(u; A), A = D t - V(u; A), 

which is conventionally called the Lax pair. Here IT, V are square matrices, whose entries are functions of 
the dependent variable u and its x-derivatives and certain rational (in some cases elliptic) functions of 
the spectral parameter A, such that equation ((3) is equivalent to the condition of commutativity of these 
operators 

[L, A] = D t (U) - D X (V) + [U, V] = 0. (4) 

The latter is often called a zero curvature representation or Lax representation of equation |(3). In this paper 
we mainly consider U and V to be 2 X 2 matrices and their entries to be polynomial in spectral parameter 
A. 

Generally speaking, symmetries of an evolutionary equation are its compatible evolutionary equations. 
Integrable equation ((3) has an infinite sequence of commuting symmetries 

u ft = T k (u,u X/ ...u x ... x ), ken, (5) 

which can be associated with a commutative algebra of linear operators 

A k = D tk - V*( U ;A), [A\A!] = Q. (6) 

Similar to equation |(3j system (0 is equivalent to [L, A k ] = 0. Operator A and equation (0 can be 
considered as members in the sequence of operators \A k ) and symmetries I0 respectively, for particular 
values of k. The commutativity of operators can be seen as a compatibility condition for the infinite 
sequence of linear problems 

0,010 = !i(u; A)W, D h (W) = V k (u; A)W, (7) 

i.e. the condition for the existence of a common fundamental solution W of all these problems, det W + 0. 
To Darboux transformations we refer as a linear map S acting on a fundamental solution 

S: ^^^ = M^, detM^O (8) 

such that the matrix function W is a fundamental solution of the linear problems 

D X (W) = U(u; A)W, D tk (W) = V k (u; A)W, (9) 
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with new "potentials" u. The matrix M is often called the Darboux matrix. The entries of the Darboux 
matrix M are a rational (elliptic) functions of the spectral parameter A. As a function of A the determinant 
of M may vanish only at a finite set points on the Riemann sphere (the parallelogram of periods). The 
Darboux matrix M depends on u, u and may also depend on some auxiliary functions g (or parameters, 
if g is a constant). For examples of such g, we refer to the matrices in (fT8b . So we denote the Darboux 
matrix by M = M(u, u, g; A). 

From the compatibility of 10 and it follows that 

D X (M) = LT(u; A)M - MU(u; A), (10) 
D tl (M) = V*(u;A)M-MV*(u;A). (11) 

Equations (fTOb and iflj} are differential equations which relate two solutions u and u of 10 and (0. In 
the literature they are also often called Backlund transformations. 

A Darboux transformation maps one compatible system (0 into another one (0. It defines a Darboux 
map S : u i-> u. The map (0 is invertible (det M + 0) and it can be iterated 

■■• = M(u / u / g;A)^^W = M(u,u,g;A)^^^ = M(U / U / g;A)W^W = M(^ / ^ / l;A)^ 7 ^ 
It suggests notations 

...vp_i = % W = =W, W 2 = VP,..., 

. . . U-i = U, Uo = U, Ui = U, U2 = U, . . . , 
• • • gl = g, gO = g, gl = g, g2 = g, 

With a vertex fc of the one dimensional lattice Z we associate variables and Ujt; with the edges 
joining vertices k and k + 1 we associate the auxiliary functions (parameters) g/t and the matrix Mjt = 
M(uk, u/t+i, git; A). In this notations the Darboux maps »S and increases and decreases the subscript 
index by one, and therefore we shall call it as *S-shift, or shift operator S. In what follows we often shall 
omit zero in the subscript index and write u, g instead of Uo, go- 

In these notations flOl and the sequence iflj} are a hierarchy of compatible systems of differential 
difference equations. When the resulting equations from fl0t and ifTTIl are in the evolutionary form they 
form an infinite dimensional Lie algebra of commuting symmetries. The existence of an infinite algebra 
of commuting symmetries is often taken as a definition of integrability of the equation (and of the whole 
hierarchy of symmetries) fT9ll20ll2Tll22l . 

In order to illustrate this construction we consider two examples: (1) the well known example of the 
nonlinear Schrodinger equation; (2) the new results on differential difference equations corresponding 
to the dihedral reduction group D2 — Z2 X Z2 (the group of Klein). In the next section we will use these 
examples for illustration of the derivation of recursion operators. 



2.1 The nonlinear Schrodinger equation 

The nonlinear Schrodinger (NLS) equation 



2p t = Pxx ~ 8p 2 q 
2q t = -q xx + 8q 2 p 



(12) 



has a zero curvature representation 10 where 



U(u;A) = I ° 2 P + A I ? , u = [l ), (13) 



2q f " \ -1 )' \ q 

The NLS equation has an infinite hierarchy of commuting symmetries. The matrix part of the corre- 
sponding linear operators A k = D tk - V k have the form 

V° = /, V k+1 = AV k + B k (u) (15) 
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where / = diag(l,-l) and B k (u) are traceless matrices with entries depending on p,q and their x- 
derivatives. Matrices B k (u) can be found recursively Il24l and in particular we have 

* = U<„;A>, V> = V<„;A) V - AV(u;A, + \{ £-_<£ ). «,6) 



Symmetries corresponding to A , A 1 , A 2 and A 3 are 

Pk = 2 P Pti = Pa- Pt 2 = |Pa.a - 4p 2 r? p, 3 = ipns, - 6pqp x 

qt = ~ 2 q qh=q x q t2 = -\q X x + 4q 2 p qt 3 = \qxxx-6pqq x . 

It is known (see 1 25 1 , 1 26 1 ) that any Darboux matrix for the NLS equation is a composition of the following 
three elementary Darboux matrices: 



(17) 



M(u, m, /; A) = I A + / \ ) , N(u, ui, ft; A) = ( ^ J J , K(u, u x ,a;A) - ( q a -i 



(18) 



and their inverses. 

The Darboux map and/or the corresponding differential difference equation for each above elementary 
Darboux matrix can be derived as follows. 

(K): The map with the Darboux matrix K does not depend on the spectral parameter A and variables 
u, ui . Moreover, it follows from l flO|l that a is a constant (does not depend on x) and p\ = a 2 p, q\ — a~ 2 q. 
Thus in this case the Darboux map is a gauge transformation corresponding to a point symmetry of the 
equation |[T2l l and its hierarchy of symmetries. 
(N): Substitution of JV(u, ui, h; A) in flO) yields 

1 

qi = -, h x = 2p 1 q 1 - 2pq, p x = -2hp, q u = 2hq x (19) 
and therefore we have got the following explicit Darboux map 

2 P(Px\ 1 

n = ?i-Z\j] x ' ^ = p 

After the change of the variables p — exp(<p) (and thus q = exp(-<p-\)) equations (fT9b yields the following 
system of evolutionary equations 

(p x = -2k, h x = 2 exp(0i - (f>) - 2 exp(cp - (20) 

which after the elimination of variable h takes form of the Toda lattice 

(pxx = 4 exp(0 — </>_i) — 4 exp((/>i - <p). 

It is an infinite chain of differential equations for dependent variables <p n , n e Z. 

Notice that <p = Inp and h = — y. Using fl9l l to eliminate x-derivatives from (Il7t one can obtain 
symmetries of the Toda chain l l20l : 

cp t0 =2 (p h =-2h cf) h =2h 2 -2(S + l)exp( ( p-(P- 1 ) 

h t0 =0 h k =2(5-1) exp(<£ -cp-t) hfc = -2(S-l)(exp(<f>-<M(fe_i + fe)) 

cp t3 = -2h 3 + 2 exp((p - (p-i)(2h + h-\) + 2 exp(0i - <p)(2h + h) 

h h = 2(S - j + h-ih + h 2 ) exp(0 - + exp(2cj> - 2^_i) + (S + 1) exp(0 - (p- 2 )) ' 
Symmetries l|2T1l have Darboux-Lax representation 

N tk -S(U*)N + NU^0 N = (J£^ ^T)' W 

and matrices U°, . . . , U 3 are obtained from V , . . . , V 3 given by |[T5l and l fT6t by elimination of x-derivatives 
using lfl9]> . In the next section we shall derive a recursion operator for generating symmetries of the Toda 
lattice. 
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(M): Substitution of M(u, ui, /; A) in ( TTOb leads to the following system of differential-difference equations: 

f x = 2p 1 q 1 -2pq, p x = 2 Pl -2fp, q x = -2q- t + 2f- X q. (23) 

Symmetries of this system can be found from the Lax Darboux representations (fTTb where V k given by 
lITBl l and (Tl6b after elimination of x-derivatives with the help of I 



Pk = 2 P Ph = 2pi - 2/p p t2 = 2(/ 2 p - / Pl - /i Pl - p 2 ^ - pp x q t + p 2 ) 

<?fo = " 2 <? ?fi = " 2i ?-i + 2 /-i<? <7i,f 2 = 2 (/-i<7 + fl ~ /V + Wi + Pi<7? " <7-i) ( 24 ) 

/fo = A = 2 p^i - 2 p<? A = 2 («s - + piq - if-i + f)pi) 

System l l23l and its symmetries ll24ll have first integral O = / - pq~\, so that O t( . = 0. Indeed, 
detM(u,Ui,/; A) = A + / - p^i should be a constant (does not depend on x,/>) because the matrices 
U, V are traceless (Abel's Theorem). Thus one can set fa - pktfk+i + <%ki where e C is a constant. We 
can eliminate / from the system 11231 and it leads to 

p. v = 2pi - 2p 2 q 1 - 2ap, q x = -2q. x + 2p_^ 2 + 2a- X q, (25) 



whose symmetries can be obtained from 1124b by the same elimination of /. When a; c = 0, equation 
becomes the Merola-Ragnisco-Tu lattice under an invertible transformation x = it, p = u and q\ = v 
listed in Section l4?10l 

2.2 Equations corresponding to the dihedral reduction group D 2 - Z 2 X Z 2 

Integrability of the equation 

2p,= p xx + iq x - 8(fq) x 

2q t = -q xx +4p x -8(q 2 p) x KM) 

has been established in [20 J. This equation can be seen as a non- trivial inhomogeneous deformation of 
the derivative nonlinear Schrodinger equation. The corresponding Lax pair 

L(u; A) = D x - V 1 {u; A), A(u; A) = D t - V z (u; A) 

has the matrix part of the form 

V\u; A) = 2pa x (A) + 2qn 2 {\) + 2a 3 (A) (27) 

V 2 (u; A) = w(A)V\xx; A) + hZJLl ^) _ £ P a2 ( A ) - 2p<7a 3 (A) (28) 

where 

JH=(S VH=^(i .!)• 

Matrices V fc (u; A) are invariant with respect to the group of transformation (the reduction group [8)) 

(5 _?)^(u;-A)(j _\) = V\n;M (° J ) V^ 1 ) ( ? J ) = VW), 

which is isomorphic to the dihedral group D2 - Z2 X Z2 (the group of Klein). Function 

A 2 + A" 2 

w{A) = (30) 

is a primitive automorphic function of the corresponding Mobius group {10(A) = w(— A) = w(A~ 1 )) and 
any rational automorphic function of this group is a rational function of zv(A) |3|. A hierarchy of higher 
symmetries of equation j26l l can be generated by the Lax operators A k - D tk — V , where the matrices V k 
are of the form 

V k+l = w{A)V k + r*ai(A) + ^a 2 (A) + r*a 3 (A) (31) 

and the coefficients r k , ti, are polynomials in p, q and their x-derivatives and they can be found recur- 
sively. 
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It has been shown in 1271 that an elementary Darboux matrix for the Lax operator ll27t can be written in 
the form 

M(p,q lr f,g;A) = f(w(A)l + a 3 (A) + pai(A) + qia. 2 (A) + gl) (32) 
where I is a unit matrix. Matrices V k are all traceless and therefore the determinant 

detM(p, qi , f, g; A) = f 2 (lw(A)(g - pq{) + 1 + g 2 - p 2 - q\) = 2ro(A)0 1 + O 2 
does not depend on x and tk- Thus we have two invariants 

& = f(g-p<n), <s> 2 = f 2 {i + g 2 -v 2 -q\)- (33) 

By choosing an appropriate scaling M — > yM we can make O 1 = 1 (if O 1 + 0) or O 2 = 1 (if O 2 # 0). 
There are three essentially different cases l27ll : 

I. O 1 = 0, O 2 = 1, in this case the determinant is a constant (does not depend on A); 

II. O 1 = 1, O 2 = +2, in this case the determinant has two double zeroes and is a square of a rational 
function of A; 

III. O 1 = 1, O 2 = 2a, the determinant has four distinct zeros in the complex plane A (assuming that 
a + ±1). 

In the case I we have 

g = m (34) 

and / can be found as a solution of equation / 2 (1 - p 2 )(l - q 2 ) = 1. Thus we can determine the matrix 
M(p, qi,f, pqi; A) in J32t . Substituting it and V x (u; A) given by l|27t instead of !i(u; A) into JlOl l leads to the 
system 

p x = 2(1 - p 2 )( qi - q), q x = 2(1 - q 2 )(p - p_ a ), (35) 
which can be written as a scalar equation for one function i>2« - Pn, Vin-i - 

i> v = (1 - z; 2 )(^ 1 - (36) 

In the case II we shall take 20 1 + O 2 = (the other choice of the sign would eventually lead to a point 
equivalent system). We have 

20 1 + 2 = / 2 (1 + g + p + qi )(l +g -p- qi ) = 

Let us choose 1 + g + p + q\ = (the second choice 1 + g — p — q\ = would lead to a point equivalent 
system). Then the substitution of M(p, q\,f, -p - q\ - 1; A) and V T (u; A) given by l|27|l instead of !i(u; A) 
into JTOt leads to the system 

p x =2(l+p)(S-l)(q-p-pq), q 1/X = 2(l + q 1 )(S-l)(p-q-pq). (37) 

In the generic case III of the elementary Darboux transformation we have 

O 2 - 2a<D 1 = (g 2 - lag - p 2 - q 2 + 2apq x + l)f 2 = 0, /= ■ 1 ■ ■ (38) 

- m 

The substitution of M(p,qi,f,g;A) and V l (u;A) given by (|27t instead of U(u;A) into (TlOb leads to the 
system 

Px = 2g(p 1 -p) + 2qi -2q- 2p(p 1 q 1 - pq), 
qi,x = 2g{qi -q)+ 2pi - 2p - 2q 1 (p 1 q 1 - pq), 
g x = 2p{p x -p) + 2qi (qi - q) - 2g{p l q l - pq), 
fx = 2f(piqi - pq). 

The invariants O 1 , 2 are the first integrals of this system. Thus the functions / and g can be eliminated 
using the first integrals 1(38)1. Note that in l l38l the parameter a is a constant in x, tk, but it may depend on 
the shift variable, so that S{a) = a\ . 
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3 Recursion operators for differential-difference equations 



In this section we show how to derive a recursion operator using a Darboux-Lax representation for a 
differential-difference equation. Our construction is a natural generalisation of the method used in the 
theory of integrable PDEs lMI2SI23l3Qj|3ll. Tne main idea of me method is based on the fact that 
matrices V k (u, A) of the operators A k = D tk - V k (u, A) corresponding to a hierarchy can be related as 

V k+1 (u, A) = ju(A) V k (u, A) + B*(u, A), (39) 

where ji(A) is a rational (elliptic in the case of the Landau-Lifshitz equation) multiplier and £>' c (u, A) is a 
rational matrix with a fixed (i.e. k independent) divisor of poles. If the system and its Lax representation 
is obtained as a result of a reduction with a reduction group G, then the multiplier /.((A) is a primitive 
automorphic function |3] of a finite reduction group, or in the elliptic case is one of the generators of 
the G-invariant subring of the coordinate ring l29l . The matrix B k (u, A) also depends on the dependent 
variables u and its x-derivatives. 

The substitution of j3"9l in the Lax representation [L, A k ] = © results in 

D tM (L) = fi(A)D tt (L) - B k L + LB k . (40) 

One can use equation | |40)| to find B in terms of variables U, U^, Uf l+1 and x derivatives of u. Then 
equation j40l can be regarded as a recursive relation Uf i+1 = ??(Uf, ), where % is a linear pseudo-differential 
recursion operator mapping a symmetry to a new symmetry. A recursion operator can be related 
to bi-Hamiltonian structure. Indeed, if Hi and H2 are two compatible Hamiltonian operators, then 
% = HiH^ 1 is a Nijenhuis recursion operator (32j [TTJ. The sufficient condition for H to be a recursion 
operator is J33J 

D tk CR) = [¥ k ,n 

where F* is the Frechet derivatives of F . 

A similar construction can be used in the differential difference case fl3l . The substitution of l l39l in the 
Darboux-Lax representation ifTTTl results in 

D fw (M) = ii(\)D tk (M) - S{B k )M + MB k . (41) 

Equation ((4TJ enables us to express the entries of the matrix B k in terms of u, Uf t , Uf t+] and their 5-shifts. It 
enables us to find a linear pseudo-difference operator % such that Uf (+1 = < R{\x tk ), i.e. a recursion operator 
for a differential-difference hierarchy of commuting symmetries. Like the differential case, if we know 
two compatible Hamiltonian operators Hi, ^ or a compatible pair of a Hamiltonian fi and a symplectic 
I operators then % = < 7Yi'7Y~ 1 and % = HI are recursion operators. From this construction it follows that 
a Darboux matrix M and a multiplier fi(A) defines a recursion operator completely and uniquely. 
In this section we shall illustrate this construction on a few examples. In the next section we shall present 
an extensive (but far not complete) list of integrable differential-difference equations with recursion 
operators, multi-Hamiltonian structure and Darboux-Lax representations. 



3.1 Differential difference equations from the NLS equation 

We are going to illustrate the construction using the Darboux matrices M and N fTHI l related to the NLS 
equation. In this case the multiplier f.i(A) = A (l5l l and the matrix B k does not depend on the spectral 
parameter A. 

(N): Let us construct a recursion operator for the Toda lattice (f20b using the Darboux matrix N l l22t and 
the multiplier ju(A) = A. We substitute 

B = (: l) (42) 



in the equation 

D tk+1 (N) = AD tk (N) - S(B)N + NB. (43) 
Linear part in A of d43l leads to the system of equations 

(S - l)a = \, b = -<fi tl exp((p), c\ = -(p tt exp(-c/)). 
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Thus we can find all coefficients a = (5 - 1) 1 h tt , b = - exp((p)(p tk , c = - exp(-<j)-i)S l (p tk of the matrix B. 
Then the A independent part of d43l l reduces to equations 



. +i = bh _ a _ ax = - {S + 1)(S - lT l h tk - h<p tk , 



K + i = ( a ~ a i) h + cexp((p) - h exp(-<p) 
which leads to a recursion 



h 



-h\ - exp(cp - cp-i)S 1 cp tk + exp(cpi - §)S<p tk , 



n 



with the pseudo difference operator 

n = 



-h -(5 + 1)(5 - l)" 1 

exp(cpi - cp)S - exp(cp - cp-^S' 1 -h 



(44) 



Starting from the seed symmetry (p to = -2, h to = one can recursively produce the hierarchy of 
symmetries j2"i~l l. 

Similar to the case of scalar discrete equations IfTTl l34l , the canonical series of the densities of local 
conservation laws can be found taking residues pu = res < R k . In the case of matrix pseudo-difference 
operator J?T a residue resJ?l is defined as follows: 

Any pseudo-difference operator J?l can be uniquely represented by its Laurent series 

00 

ji=Y J w m - k s'"- k , 

k=0 

then the residue is defined as resjft = trace(J?T ). For example, we rewrite the recursion operator (f44b as 



It follows that 











exp(0i - <p) 



5 + 



-h 



-2 
■ exp((j) - (p-{) 



s- 1 



pi = res?? = -2h, p 2 = res?? 2 =2h 2 -2 exp{$i -(p)-2 exp((p - cp-i), . . . , 

are conserved densities for the Toda lattice l l20l l. Indeed, we have 

D xPl = -4(5 -l)exp(<f> -c/)-i); 

D x p 2 = 4(5 - 1) (exp(cp - + hj) . 



(M): Now we take the Darboux matrix M JT~8l > and the matrix B of the form 
A of the equation 

D tM (M) = AD tk (M) - S(B)M + MB 
it follows that 

a = (5- 1) _1 A, b = ~Ph, c = q tk . 
The A independent part of d45b leads to the recursion relation for the system ( f23t 



From the linear part in 
(45) 



( ft M ' 




( -fK+S(qPh) + vqt k , > 




{ k ) 






-2 V {s-iy l f tk -f tkV -f Vtk+ S{ Vtk ), 


= n 


Vk 






{ 2q{S - - qS^iftJ + 5^ - f-iq h j 







(46) 



where a pseudo-difference recursion operator is of the form 





{ -f 


q\<S 


V > 






\ 


n = 


-p 


-f + S 





+ 




-2p 




l -qS- 1 





-A+5- 1 , 




I 


2q j 



(5 - 1) -1 (1 0). 
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Using the first integral of the system (f23l l and eliminating fa = pktfk+i + &k we obtain the corresponding 
recursion relation and operator for the system d25t as follows: 



Pt k+ 1 \ _ CQI I Vh 



K'\ ™ , (47) 



.S-lpqt-a -fS-2pp- 1 \ I -2p , t v 



where 



As we point out that equation I l25ll is related to the Merola-Ragnisco-Tu lattice, we can check that we 
obtain the same operator from the recursion operator for the Merola-Ragnisco-Tu lattice in Section 14.101 
using formula 



3.2 Difference equation corresponding to the dihedral reduction group 

In this section, we demonstrate how to compute a recursion operator for equation l l37t , whose Lax 

representation is invariant under the dihedral reduction group. 

After a simple change of variables p— >p — 1,0— »0 — 1 equation l|37t becomes 

p x = 2p(S-l)(2q-pq), q u =2q 1 (S-l)(2p-pq). (48) 

This equation is the relativistic Volterra lattice in Section l4~9l under the scaling transformation 

t 

p = -2v, q = -2u, x---. (49) 
The corresponding Darboux matrix and the multiplier are 

A 2 + l-p-fl A(p - 1) + A" 1 ^! - 1) \ „ m _ l n2 ,-2, 



M = f\ it i \ i-l, ^^ i-2 i / M(A) = ^(A z + A" z ). (50) 

This information is sufficient to find a recursion operator, equation (T48b itself and its hierarchy of local 
symmetries. It follows from 

detM = -(A - A" 1 ) 2 / 2 P'?i 

that f 2 pqi does not depend on x (Abel's Theorem). Thus we can set / = {pq\)~^- . 
Expression l|3il suggests that 

B = a&i + b&2 + ca3 . 

Now in equation (f4lTl the right and left hand sides are rational matrix functions in A. The left hand side 
has only simple poles. Vanishing the coefficients at the third order poles in A -3 in the right hand side is 
equivalent to 

(q 1 +l)pSq fk +(l-q 1 )qip, l! \ D1 > 



-- 2 t(c + ci)(p-l) + 
h = l(c + c 1 )(q 1 -l) + ~ 



The second order poles vanish iff 



c ~ Cl = — ^ — • ( 52) 



Now we can simplify d5Tb by eliminating c\ and c_i from a, 

= c(p - 1) - |p ft , 

b = c((7 _ 1)+ j (7ft . ^ 
The first order poles in equation dill vanishes if 

Piiw = C P(S - - 2«) + (p - 0! - 1 - \p\q\)p h - \pq x Sp tk + 

(p - \p 2 )qh + (p - 2-ppi)Sqt k , 

(54) 

<7fw = c <?(3 - l)(P-i«-i " 2p-i) + (0 + p-i - 1 - \p-iq-{)q k - \p-iqS l q tk + 

(q - \q 2 )Ph + iq - \qq-i)S~ l Ph- 



li 



The A independent part of equation d43~b is satisfied after the substitution of d52l l, d53l l and ( f54t in <f4T)) . 
Equation d54b together with d52l l is a recurrence relation. The differential-difference equation (a seed 
symmetry) p x = 2K^\ q x = 2 K^ ll48l l can be recovered from this recursion by taking the derivative of the 
right-hand side of equations ll54t with respect to c, that is, 

<$\ K ( V) = (p(S - l)(pq - 2q), q{S - l)(p-xq-x - 2p_i))- 

It is not surprising. Indeed, in order to solve equation ((52) to find c we need to invert the difference 
operator S — 1, whose kernel is a field of constants. As a result, the vector p x , q x may contribute to 
P't+i ' Ih+i w ith an arbitrary constant coefficient. 

We now write explicitly down the recursion operator corresponding to the recursion ((54) as follows: 



n 



p + qi-l- \pxq x - \pqiS p-\f + {p- \ppx)S 

q-\q 2 + {q- \qq-\)S~ l q + p-x-1- \p-\q~\ - \p-\qS~ x 

I p(S-l)(pq-2q) ) (5 _ ir i ( l J_ 5). (55 ) 
\ q(S - l)(p-xq-x - 2p- x ) J v 2p 2q x 

Obviously K + I, where I is the 2x2 identity matrix, is also a recursion operator, which can be recovered 
from the recursion operator for the relativistic Volterra lattice listed in Section l4~9l using the transformation 



4 A List of Integrable Differential-Difference Equations 

In this section, we present a long list of integrable differential difference equations with their Hamiltonian 
structures, recursion operators, nontrivial generalised symmetries and their Lax representations. To be 
self-contained, we introduce notations and remind some definitions of the objects in our list in terms of 
the Lie derivatives. The theoretic background and detail definitions of the Hamiltonian and symplectic 
operators can be found in monographs l35l l36l . Here we would like to mention that there are some 
recent developments in the theory of nonlocal Hamiltonian structures for nonlinear partial differential 
equations in J37). 

Let u = (m 1 , . . . , U ) be a vector-valued function of variable neZ and time variable t. An evolutionary 
differential-difference equation of dependent variable u is of the form 

u f = K[u], (56) 

where K[u] means that the smooth vector-valued function K depends on u and its shifts u, = S'u, i € Z. 
In all our examples after an appropriate point change of variables we can consider K[u] € T~ N where 
f = (C, u, S) is a difference field of rational functions of {u\ \ k e Z, i — 1, . . ., N}. A Frechet derivative a* 
of a e f is defined as a row vector of the difference operators 

L. da da „ k 

M < } - 

feZ k k 

thus the Frechet derivative of K[u] is a difference operator with square matrix coefficients and the entries 

of the coefficient matrices are elements of f. 

A variational derivative of a e f is a column vector 

M4 -"'.(I) = (|r ^)2>' ( * 

v ' feZ 

If equation d56l l is Hamiltonian, then it can be written in the form 

where "H denotes a Hamiltonian (pseudo-)difference operator and / is a Hamiltonian function (or the 
Hamiltonian of the system). 
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Definition 1 Given a differential-difference evolutionary equation (56$ , we define that 

• G is its symmetry z/LrG := [K, G] := G*(K) - K*(G) = 0; 

• 'H is its Hamiltonian operator if a Hamiltonian operator *H satisfies 

Lk'H := W*[K] - K*<W - <H< = 0; 

• J is zfs symplectic operator if a symplectic operator I satisfies 

LkT := J*[K] + K + + J + JK* = 0; 

• <R is its recursion operator ifL K ft := 7?*[K] - K*ft + 7?K* = 0. 

Here Lk denotes the Lie derivative, ★ means the Frechet derivative and t means the formal conjugation of the 
difference operator. 

If a symmetry of equation d56t is explicitly dependent on u, with i + 0, we say it is a generalised symmetry 
An equation is integrable if it possesses infinitely many generalised symmetries depending on finite sets 
of variables u„ whose sizes are increasing. Symmetries of an integrable systems can be generated by a 
recursion operator which is often nonlocal. Sufficient conditions on nonlocal pseudo-difference recursion 
operators which guarantee the production of an infinite hierarchy of commuting local symmetries have 
been discussed in IITTll . In the list, we also give partial results on master symmetries (see I38U39I for more 
details on master symmetries). 

We now look at how the recursion, Hamiltonian and symplectic operators change under transforma- 
tions (difference substitutions). If an evolutionary difference equation J56l l is connected by a difference 
substitution u = F[v], F[v] e f with another equation of the form 

v f = G[v] (57) 

then recursion, Hamiltonian and symplectic operators ("R, *H and I) for the equation | |57|| can be expressed 
in terms of the corresponding operators {K, "H and I) of the equation l l56t 

^ = F; 1 o??| u=F[v] oF* / 'H = F; 1 o'H| u=Ftv] oFr 1 , J = FtoJ| u=F[v] oF* (58) 
where o denotes the composition of the operators. 

Example 1 It is known that the Volterra equation u t = u{u\ - U-\) (section WD can be related to the modified 
Volterra equation 

v t = v\ Vl - v-i) (59) 

by a difference substitution (a Miura type transformation) u = F[v] = vv\. Operator *Hj = u{S - S~ x )u is a 
Hamiltonian operator for the Volterra equation. Notice that the Frechet derivative for vv\ is F* = Vi + vS, and 
we have F^u = v(l + *S) _1 and uF\ 1 = S(l + *S) -1 c. Now we can find the Hamiltonian operator < H\ for the 
modified Volterra equation (59$ 

% = (di + vsy l w l (s - s~ l ) w-i{vi + vsy 1 

= z;(l+5)- 1 (5 2 -l)(5 + l)- 1 u 
= v(S - 1){S + l^v 

In the same way, we can compute the second Hamiltonian operator e H% and the recursion operator %for the modified 
Volterra equation {compare with section WlS . 

In Section|2]we discussed that the compatibility of a Darboux map ||8]l with the Lax operator 10 

5(d)) = MO, Df(®) = !i(u;/\)<5 (60) 

yields O 

D t (M)=S(U)M-MU, (61) 



13 



which is equivalent to an integrable system of differential-difference equations. In literature the com- 
patibility condition 1I6TI1 is often called a zero curvature representation or Lax representation of equation ll56l l. 
Since it involves a Darboux matrix and a Lax operator, it is more appropriate to call it a Darboux-Lax 
representation. However, we still call it a Lax representation to be consistent with the literature. In the 
following list, we simply give the expressions of both matrices M and U for Lax representations. 
There are many publications where an integrable system emerges as a compatibility condition of two 
linear problems with scalar linear difference operators L and A 

L(p = A(p, <p t = A<f>, (62) 

where <p is an eigenfunction of L corresponding to eigenvalue A and At = 0. Equation d56") is equivalent 
to their compatibility condition 

D f (L) = [A, L]=AL-LA. (63) 

This approach really resembles the original Lax formulation where differential operators are just replaced 
by difference ones. In the theory of ordinary differential equations a scalar high order equation can be 
represented in the form of a system of first order equations. Similar we can do in the case of scalar 
difference equations of high order and represent them by a first order difference systems. In this 
way we can re-write a scalar representation (62), ((63) in the form of a first order matrix Darboux-Lax 
representation ((60), ( |6Ti . 

Example 2 For the Volterra Chain listed in section \4~l\ its matrix Lax representation is obtained from the scalar 

Lax representation with L = S + uS~^ and A = S 2 + u\ + u. 

Let <I> = (cp 1 , cp 2 ) T = {cb, -5 -1 (/>) T . We can rewrite Lcp = Acp in ([62]) as 



S(0) = S 



S<M_/ A u\( ch 



that is, Sty 1 ) = A(p x + u<p 2 and S(<p 2 ) = -cp 1 . Now cp t = Acp can be rewritten as 

(p) = S 2 (p : + (wj + u)^ = SiAcp 1 + ucp 2 ) + (mi + m)^ 1 = AScp 1 + ucp 1 = (A 2 + u)<p l + Aucp 2 



and 
Thus 



1 ^ ' A 2 + u Au w Al 



In our list we present either scalar or matrix Lax representation. 
The reader should bear in mind that 

• Relation between nonlocal difference operators, e.g. l l67t , should be understood as identities in the 
non-commutative field of pseudo-difference Laurent series. 

• In the computations with pseudo-difference operators, we often use identities similar to "integra- 
tion by parts": 

(S - irVi - /)(«S - l)" 1 = f(S - I)" 1 - (S - I)" 1 /! 

• Since for any constant c e C we have (S - l)c = 0, the action of operators involving (S - l) -1 
is not uniquely defined. The corresponding results given in the list are up to these "integration 
constants". 
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4.1 The Volterra Chain 

• Equation IflOl : 

u t = u(u\ - u-i) (64) 



• Hamiltonian structure 14011411 : u t = Hi5„f 

"Hi = u(S - S'^u, fx = u 



9^2 = u(SuS + uS + Su - nS~ Y - S~ Y u - S~HS~ )u, f 2 = — In u 



Recursion operator: 



<R = <H 2 <H- 1 = uS + u + ui + uS' 1 + u{u x - U-\){S - l)" 1 ^ 



= u(S - S^)u (-(S - l)" 1 + S(S - l)" 1 -) 
\u u I 

• Non-trivial symmetry IHT1I4211 : 

9?(Mf) = U [u\U 2 + U 2 + UU\ - UU-\ - U 2 _ x - U-\U-2^ 

• Master symmetry I4lll38ll43l: 

Kill) = nu t + u(2u\ + u + U-i) 

Lax representation | 



L = S + uS' 1 , A = S 2 + mi + u 
It can also be written in the matrix form 



A u \ i j I A 2 + u Aw 



M i-i o y u =\ -a »_ a j < 65 > 

This equation is also known as Lotka- Volterra model, the Kac-van Moerbeke lattice or the Langmuir 
lattice l45l . The so-called Kac-Moerbeke-Langmuir equation t40l 

zv r = w(p\ - w e _^), e + is a constant. 

is related to l l64t by the point transformation u = vf and t = ex. This equation is also written as 

Wt = exp(w + W\) - exp(zv + W-i), 

which can be transformed into l l64l l by the transformation u = exp(it> + u>\). 

4.2 Modified Volterra equation 

• Equation f46l[42l: 

U t = U 2 {U\ - U-i) 

• Hamiltonian structure l47ll42l : 

■Hi = u{S - 1){S + iy l u, fx = uui 
f H 2 = u 2 {S-S- l )u 2 , f 2 = hxu 
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Recursion operator: 

<R = <H 2 , H; 1 = u 2 S + 2uui + m 2 5 _1 + 2w 2 (ui - u-i)(S - l)" 1 - 

u 

Non- trivial symmetry Il42l : 

K{u t ) = U 2 u\(ll2 + U) - l^U^U + M_ 2 ) 



Master symmetry 



• Lax representation 



<R{-) = nu t + y(3ui + w-i) 



u | 1T / uu-i Au-i 
M = l , ; U = \ , l2 

-u A I \ -Au A + uu-i 



The Modified Volterra equation is also known as discrete modified Korteweg-de Vries equation. Under 
the Miura transformation w = UU\ it can be transformed into the Volterra Chain w t = w(w\ - W-\) as in 
section |4~T1 



4.3 Yamilov's discretisation of the Krichever-Novikov equation 

• Equation 1351: 

^ _ R{ui,u,u-i) 

Ml - M_i 

where R is the polynomial with constant coefficients a, jS, y, 5, e defined by 

R(u, v, w) = (av 2 + 2f3v + y)uw + (/3c 2 + Av + 5)(u + w) + yv 2 + 25v + e . (66) 



Two non-trivial symmetries t49ll50ll5lll34l[TTll52l: 

- U-2 I 



(2) _ R(u, u-i, u)R{u lr u, «i) / 1 



(U\ - U-i) 2 \«2 - U U 



R{3) = R(mi,m,mi)R(m,M-i,m)/ | ^W\ +K(1 ) K (2)/ 1 , 1 \ 

(«i - M_i) 2 \(«2 - U) 2 (U - U-2) 2 ) \U2~U U-U-2/ 



Hamiltonian structure IHT1I52I : 



where 



A 
A 



<H=AS-S~ 1 A+2 K m (S - l)" 1 ^ 2 ' + 2 R (2) (S - l)" 1 ^; 
<H = AS 2 - S~ 2 A + BS-S~ 1 B + K {2) (S - iy\S + 1)K {2) 
+ 2 K m (S - ly^SK® + 2 R (3) (S - 1)- : X (1) , 



R(m 2/ U\, u 2 )R(ui, u, u\)R(u, u-i, u) 
(ui - u_i) 2 (w 2 - u) 2 

R(ll3, U%, U^)R{ll2, Ml, « 2 )R(Mi, M, U\)R(u, U-i, u) 



(Mi - M_i) 2 (M 2 - m) 2 (m 3 - Mi) 2 

B = 2a{ K(1) - duR ( Ul ' u ' u ^ + d z R{ux,u,ux) \ + 2R(m / M-i,m) ^, Mi) K ( 2 )\ 

' M - M_2 2(Mi - U-\) Adudui ) (Mi - M_i) 2 ^ ' 
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Symplectic operator: 



I 



1 



R{ut,u,ui) 



R(ui,u,Ui) 



• Recursion operator IITT1I52I : 



R = <HI and H = <HI 



The recursion operators H and H satisfy the algebraic equation 

(2fc - hf = 4(7? + hf - g 2 (R + 1 2 ) - g 5 , 



(67) 



where h,h,g2 and g3 are the relative and modular invariants related to h = R(u\,u,u\) and a quartic 
polynomial f(u) = (d Ul h) 2 - 2hd 2 li h defined by 

S2 = ^(2ff IV -2f'f"'+(f") 2 ), 

= ^ (l2//'7 TV - %f'ff V - 6/(/"') 2 + 6/7"/'" - 2(/") 3 ) , 



^3 



h = I (hd 2 u d 2 u h - (d u h)(d u dlh) - (d Ul h)(dld Ul h) + (dlh)(dlh)) + ^{d u d lh hf, 



h = - det 



( h d u h d\h 

d u Ji d Ul d u h d Ul dlh 
[dlh d 2 Ul d u h dldlh 



4.4 Integrable Volterra type equations 

The classification of integrable Volterra type equations of the form 

u t = f(U-i,U,Ui), 

where / is a smooth function of all its variables was obtained by Yamilov using the symmetry approach. 
In his remarkable review paper [42J, he presented the following complete list of integrable Volterra-type 
equations (with higher order conservation laws) up to the point transformations: 



VI. 
V2. 

V3. 

V4. 
V5. 
V6. 

V7. 

V8. 

V9. 

V10. 

Vll. 



Ut = P(ll)(U\ - U-i) 



III 



= iV) p !_) 

\U\+U U + U_i / 

-!-) 

- U-i) 



u t = Q(u) 
"I 



III - u u 

R(ui, u, U-i) + vR{u\, u, u\) ll2 R{u-\, u, U-\) 112 



U\ - U-i 

u t = y{u\ - u) + y(u - u-i), y' = P(y) 

P(y) 



ve {0, ±1} 



u t = y{u\ - u) y(u - U-\) + [i, 



y 



(U e C 



iii 

u t = 
u t = 
u t = 



y(u\ - u) + y(u - U-i) 

1 

y{ii\ + u) — y(u + u-x)' 
y(ii\ + u) - y(u + ii-i) 
y{ii\ + u) + y(u + U-i)' 
y{u\ + u) + y(u + m_i) 



y'=P(y 2 ) peC 



y' = Q(y) 
P(y 2 ) 



y 



y(u\ + u) - y(u + U-\)' 
(1 - y(«i - »))(! - y(u - M-i)) 
y(u\ -u) + y(u - u-\) 



y 

Q(y) 
y 



P(y 2 ) 



, neC, 



(68) 
(69) 

(70) 

(71) 
(72) 
(73) 

(74) 

(75) 

(76) 

(77) 

(78) 
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where P and Q are polynomials with constant coefficients a, j8, y, 5, e defined by 

P(u) = au 2 + fin + y, (79) 
Q(m) = «m 4 + y3w 3 + yu z + 5u + e, (80) 

and the polynomial R is defined by ||66}. As stated in the paper l42ll , the problem of constructing 
the generalized symmetries for all equations (Vl)-(Vll) remains open although the master symmetries 
for some forms of equations in the list are known till l38l . We know that the Miura transformation 
u = y(u\ - u) transforms equations (V5) and (V6) to (VI) and equations (V7) and (Vll) into (V2), and 
the Miura transformation u = y(u\ + u) transforms equation (V9) to (V2) and equations (V8) and (V10) 
into (V3) l42l . In the follows, we present the recursion operators, Hamiltonian operators and master 
symmetries for the first four equations. The corresponding operators for other equations can be obtained 
via the Miura transformations. 

4.4.1 VI equation GOD 

Hamiltonian structure: 



<H = P(u)(S - S _1 )P(«) 

Symplectic operator: 

a(S - S- 1 ) + (au, +f} + au-i)S(S - l) -1 ^ + ^(5- irV"i + /S + ow-i) 

l (U) l (U) 

Recursion operator: 

<R = P(u)S + laum + p(u + mi) + P(w)vS- 1 + u t (S - 1) _1 ^t 

l (u) 

Non-trivial symmetry : 

Hiiit) = P(u) {P(ui)u2 + auu\ + fi(u + - P(m_i)u_ 2 - auu 2 _ t - jS(w + U-\)u-\\ 
Master symmetry: 

nu t + P(u){ai\ + — + (2 - c)w_i), ceC, when a + 0; 
a 

nu t + P(u)(cu\ + u + (3 - c)u-\), c eC, when a = 

Lax representation: 

- The case a = ji = is a linear equation. 

- The case a = 0, f> + reduces to the Volterra equation d64b by linear substitution u i-» ^~ l (u-y). 
Thus it has the Lax representation (T63b . 

- In the case a + the following linear substitution t i-> a~ l t and u i-» m — j3/2 transforms the 
polynomial P(m) in VI to the form P(w) = m 2 + c, where c = y/a - |S 2 /(4a 2 ). The corresponding 
Lax representation is of the form Il25l 

. cA _1 u \ u - ( c2 ^~ 2 + uu -\ cA~ l u + Am_i 
-u A J' \ -cA _1 m_i - Au A 2 + uu-\ 



This equation includes both Volterra chain in section |4T1 and modified Volterra equation in section [4721 
For nonlinear equation, that is, af> + 0, its Hamiltonian / depends on the coefficients in polynomial P 
defined by ((791 . When a + 0, we take f = lnP(w), otherwise, we take / = |. 

Its Hamiltonian operator, symplectic operator and recursion operator satisfy the following unexpected 
relation 

HI = aK 2 + fiK + 2y(|S 2 - lay). 
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4.4.2 V2 equation <69) 

• Hamiltonian structure: 



U\+U U + U-i U + U-i U\ + M 



Symplectic operator: 



(tt + Ml) 2 (M + Mi) z 



-(^ 2 -4ay)^r(5-l)(5 + l)- 1 r 



where p — j In ^prrj an d thus 



P(u 2 ) v 7 P(u 2 )' 

, ann fnnc 

P(«2) 



1 1 P(u 2 )' 

5 u p = + P(u 2 )' = 2u(2au 2 + B) 

K ui + u u + u-i 2P(w 2 ) v ; v y ' 



Recursion operator: 



(«i + u) 2 \(w + «i) 2 (m + mi)(u + m 1 ) (w + «_i) 2 

+P(" 2 )' (f - — ) + , P( " 2 \ 2 ^" 1 + 2u t (S - l)-%p 
\2u u + ui) (u + u-i) 2 r 

Non-trivial symmetry : 

P(u 2 )P(« 2 ) i i i \ P(w 2 )P(w 2 1 ) 



l i \ fW(«y/ i _ i \ 

+ W2 U\ — U-\) (ll + W-l) 2 VM-i+M-2 U-\—U\) 



(U + U\) 2 \U\ 

auhi 2 ! + Bu\U-\ + y 

+2P(u 2 ) 1 -1 - - 

(u + — u~i)(u + U_i) 



Master symmetry: 



P(w 2 ) 

WW* H aw 3 - jSw, when y = 0, 

U + M_i 
P(« 2 ) /3 

mm* H au 3 - — w, when B 2 - 4ay = 0, 

u + u-t 2 K r 

P(u 2 ) 

nu t + , when a = 0. 

u + M-i 

Here we only found the master symmetries for some special cases. For the given Hamiltonian operator 
above, we have 'Hb u p = < R{ut). Besides, the Hamiltonian operator, symplectic operator and the recursion 
operator satisfy the following relation: 

<H I = n 2 - 2B K + ft - 4ay. 

The Calogero-Degasperis Lattice 11381 

u t = -(l-u 2 ){b 2 -a 2 u 2 ){— — ) 

4 ' ui+u u + u-i 

is one of special cases of the V2 equation. The authors of [38] presented its master symmetry by 
introducing the time dependence for the coefficients a and b, which is in different form from what we 
did. 
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4.4.3 V3 equation (70) 

• Hamiltonian structure: 

a, Q(U) c Q(«) Q(u) Q(m) ! 

71 = o S + u t (S + IMS - 1) Wf 

U\-U 11- U-i u - U-i U\ - u 

• Symplectic operator: 

I = , 1 , 2 S - S~ \ 1 - 5„p(5 + 1)(5 - l)- 1 ^ 
(wi - u) 2 (Wi - U) z 

2 2 

-(lay - p 2 /A)-^—(S - 1)(S + l)" 1 ^ - ^t^t(«S - 1)(5 + l)" 1 



where p = \ In 77^2 and thus 



Q(«) v Q(«) r Q(u) 1 Q(u) 

+(2a6+M( ^_ (S+ir .^__^ (S+ 

anj'l rhiu 

(«i-») 



Sn 1 1 4. Q '( U) 

U\ — u 11- U-i 2(3(«) 

• Recursion operator: 



(ii\-u) 2 \(u-U-\) 2 («i-w)(m-w 1 ) (ui-u) : 

Q'(u) _ 2 , c-l o /c t\-U 

law" - Bu -y + —S - 2u t (S - 1) L 6 u p 

U\ — u (u- U-iY 



Non-trivial symmetry : 



v Q(«)Q(m) ^ Q(m)Q(«-i) 

* (M,) = (» 1 -»)2( M2 -« 1 ) + (u-,-:) 2 ^-!-^) + aQ( " )(Ul - U - l) 

+Q(«) ( 7 — ^ — ? + 2«» 2 + (— + 

- U)(U - U-i) )\U\-U 11- U-i ) 

• Master symmetry: 

Q(u) , 7 
nu* 1- au + Bu + yu, when = e = 0, 

U + U-i 

Q(u) , a _ 

nw t , when a = B = 

M + U-i 

Similar to V2 equation in section l4.4.21 we have not found the master symmetry valid for the polynomial 
Q (|80l l with arbitrary coefficients. However, if we look for the master symmetry of the form 

Q(«) 



ntt f h 7 c,w , a 6 C. 

!=0 



we can determine the constants c, for certain polynomial Q. Two examples are listed above. 

Notice that 'Hbup = ??(Uf) and the product of the Hamiltonian operator and symplectic operator gives 

rise of the square of the recursion operator, that is, 

<hi = n 2 . 
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4.5 The Narita-Itoh-Bogoyavlensky lattice 

• Equation 



U t = nC£j ii k - ^ "-*)/ V e N 

k=\ k=l 

• Hamiltonian structure llT5l : 

v v 
<H = uij^ S i - ^ S~ l )u, f = u 

i=l !=1 

• Recursion operator 1TT31 : 

^P 

<R = u(S - S-P)(S - l)" 1 [|(vSP +1 -'m - uS-'^Sr-hi - uS'T 1 , 

i=i 

where the notation FJ ._? is denoted the order of the value i, from 1 to p, that is, 

• Non-trivial symmetry: 

K{u t ) = m(1 - v S" ( p +1) )vS 1 -p Y" MyM, +p 

0<!</'<2p-l 

• Master symmetry fl3l : *R(m) 

• Lax representation l44ll : 

L = S + nS-r, A = (L ( p +1) )> 0/ 

where > means taking the terms with non-negative power of S in L' p+1 ^ . 

For p = 1,2 or 3, a few higher order symmetries are explicitly given in ITT51 , where the authors also 
studied their Hamiltonian operator, recursion operator and master symmetry for p = 1,2. 
The Narita-Itoh-Bogoyavlensky lattice is known as an integrable discretisation for the Korteweg-de Vries 
equation. It can also be presented as 

P P 

V t = v(Y[Vk ~Y\ v -k) 

k=l k=l 

which is related to the Narita-Itoh-Bogoyavlensky lattice via the transformation u = Yl^ = o v k for fixed p. 
Taking p-1, we get the well-known Volterra chain in section 14.11 Thus they can be regarded as the 
generalisation of the Volterra chain. 

Let u = Yl P k=0 u>k- Then w satisfies the so-called the modified Bogoyavlensky chain 

v v 

W t = ™ 2 (Y\ W k - \\ w -k)- 
k=\ k=l 

The recursion operator given above for the Narita-Itoh-Bogoyavlensky lattice is highly nonlocal (so is 
the master symmetry). Recently, Svinin [55] derived the explicit formulas for its generalised symmetries 
in terms of a family of homogeneous difference polynomials. The properties of these homogeneous 
difference polynomials |56| enable us to prove the locality of its infinitely many symmetries f!3l . 
A family of integrable lattice hierarchies associated with fractional Lax operators was introduced by 
Adler and Postnikov l57ll58l . One simple example is 



p-i p-i 

9 /i r l r 

u t 

k=l k=l k=l k=\ 



= u 2 (Y\ Uk - Yl u-k) - u(Y\ u k - Yl u-k), 2 < p e N, (81) 



21 



which is an integrable discretisation for the Sawada-Kotera equation. It can be considered as inhomo- 
geneous generalisation of the Bogoyavlensky type lattices. The problem to construct the Hamiltonian 
structure and recursion operator for such family of equations is still open. 

Even for in the scalar case, the classification of higher order of integrable evolutionary differential- 
difference equations is still open. 



4.6 The Toda Lattice 

• Equation l59l : 

q tt = exp((j! - q) - exp(^ - q-i) 

In the Manakov-Flaschka coordinates |60, 45 1 defined by u — exp(q\ - q), v - q t/ it can be rewritten 
as two-component evolution system: 



Hamiltonian structure 



Recursion operator: 



Uf - u(v\ - v) 

V t = U - U-i 



( U(S-1)\ r_ 1l+ V 2 

_/ uiS-S-^u u(S-l)v \ 

v(l-S-')u uS-S-hi ' fl ~ v 



(82) 



Non-trivial symmetry: 

'R 



u t \ _ I u(v* - V 2 + III - U-i) 
V t J ~~ I ll(V\ +V) - U-\(V-i + v) 



Master symmetry 



u \ ( nu t + \uv\ + \uv 

| / \ HVt + U + U-i + y 



• Lax representation: 



The Hirota nonlinear Lumped Network equation 16T1I 

U t =V\—V 
V t = V(ll - U-i) 

is related to the Toda lattice ll82l by a simple invertible transformation. Namely, let u = q and v = p_j. 
Then the variables p and q satisfy the Toda equation. All its properties can be obtained via those for the 
Toda lattice. 
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4.7 A Relativistic Toda system 

• Equation 11621 : 

cjtt = qtq-ir. 



exp(<?_i - q) 



qtq\r. 



exp(q - q t ) 



1+expiq^-q) l+exp^-^) 
Let us introduce the following dependent variables l63ll64l 

q t exp{q-qi) q t 



u = 



1 +exp(q-qi) 
Then the equation can be written as 



v = 



1 +exp(^-^)" 



U t = ll(U-i - U\ + V - V\) 

v t = v(u-i - u) 



(83) 



• Hamiltonian structure 







^2 = 

• Recursion operator 



u(l - S) 
iS-S~ 

u(l -, 

viS' 1 - l)u 



= ( (tS -i _ 1)M uS-S-'u)' A = 2 iu +v) + UV + UlU + UVl 
uiS' 1 - S)u u(l - S)v 



fl = U + V 



uS + u + v\ + ii\ + nS 1 - u(v - v\ + u-i - ii\)(S - 1) 1 ~ uS + u 
v + vS' 1 -v(u-i -u){S-l)- l \ v 

uS + u + v\ + ii\ + uS~ l uS + u 
v + z\S _1 v 



(S-l)-'( \ 0) 



Non-trivial symmetry : 



UU-i(u + U-i + U-2 + 2v + V-\) - UU\{U2 + U\ + U + 7.V\ + V2) + u 2 (v - V\) + u(v 2 - C 2 )\ 
VU-i (U-2 + U-l + V + V-\) - uv{u\ + U + V\ + v) J 



Master symmetry 



n 



• Lax representation: 



M = 



—nut + u(v + 2v\ + u + 1u\ + u-i) 
-nvt + v(n + v + u-i) 



), ((S-l)-h=n) 



Av — A 1 u-\ 
-1 



U = 



-A 2 - u-\ A l u-\ 



-U-2 - V-i 



As noticed in |64], the inverse of this recursion operator H is also weakly nonlocal: 



<r- 1 = <n x <n- 



X _*5+«_2« 
-5-U-1 45 + 5-1^+2^ + 1 



V V-i Vl I 

However, recursion operators K and ( Rr l have different starting point, i.e., seeds. For %, it starts with 

/ --" \ 

acting on the right hand of the equation while the seed for "K is a = Z\ _ v t , ■ Moreover, K acting 
on a and fC 1 acting on the right hand of the equation do not give rise to new symmetries. 
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There is a Miura transformation u - --r-h,V = — 4- between the flow corresponding to a and the 
equation, where u' and i>' denote dependent variables for a. 

In [63], the author studied other integrable equations related the relativistic Toda lattice. For example, 
the equation 

qtt = (J-ii exp(q_i - q) - exp(2q_i - 2q) - q u exp(q - q t ) + exp{2q - 2^) 
can also be rewritten as system d83l l by setting 

u = exp(q - qi), v = q t - exp(^_i - q) - exp(q - q{). 

4.8 Two-component Volterra lattice 

• Equation EC]: 



Mi = ll(V\ - V) 
v t = v(u - W_i) 



• Hamiltonian structure [40 J : 



nj \ u(S-l)v\ , 

^=[v(l-S-^u J' f' = U + V 

fj, _( u(Sv - i\S _1 )« u(uS -u + Sv — v)v \ _ 
Hl ~\v(u- S-hi + v- vS-^u v(uS - S-^v )' h ~ ln " 

Recursion operator: 

u + v x + u{v x -v){S- l)" 1 1 uS+^+ u{vi - v){S - l)" 1 i 



v + vS~ l + v(u - u_!)(vS - l)" 1 ^ u + v + v(u — u-i)(S - 

u + vi uS + ^\ I u(vi-v) L_ irl / i i\ 
v + vS' 1 u + v )^~{ v(u-u-i) j v ° ' V» 

Non-trivial symmetry: 



<R> Ut 



U 2 {V\ - V) + U (v\ 2 - V 2 + V\U\ - VU-ij 
Vt ) \ V 1 {u - U-i) + V{ll 2 - U 2 _ 1 + UV\ - U-\V-\) 



Master symmetry 



Lax representation 



u \ _ I 2nu t + u 2 + 3uv\ 
v J \ 2nv t + vu-\ + 2uv + v 2 



(84) 



L = AS 1 + v + u-i + A 1 uvS, A = A 1 uvS 

This system comes from the Volterra chain in section |4~T1 written in the variable w, that is, 

w t = w(wi - w-i), (85) 

by renaming u(n, t) = w(2n, f) and v(n, t) = w(2n - 1, f). It is related to the Toda equation 1 (82) , written in 
variables u and v, by the Miura transformation l66l 

u = uv and v = + v. (86) 

In fact, the (master) symmetries, conservation laws and local Hamiltonian structures of this system 
can be easily obtained from the Volterra chain in the same way. For instance, we can derive the first 
Hamiltonian operator 'Hi as follows: 
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A symmetry flow of the Volterra chain |85) is 

w T = w(S - S~ 1 )wQ[n] = ww\Q[n + 1] - ww-\Q[n - 1], 

where Q[n] is the variational derivative of a conserved density for l l85ll . We now write down both even 
and odd chains and rename them for the variables u and V accordingly. We have 



u T = uv\Q[2n + 1] - uvQ[2n - 1] and v T = vuQ[2n] - vu-\Q{2n - 2], 



that is, 



u(S-l)v \ Q[2n] 

vil-S-^u II Q[2n-1] 



Using the same method, we can derive Hz in the list. 

Such construction is valid for all scalar equations. It is tricky if the operator is nonlocal. 



4.9 The Relativistic Volterra Lattice 

• Equation I65ll66ll27ll: 



U t = U(V - V-\ + UV - U-\V-i) 
V t = v(ll\ — U + U\Vi — uv) 



Hamiltonian structure 





v(S - l)u 



m(1 - S-^v 




fi = u + v + uv 



(j^ _ I uv{l + u)Su-uS 1 uv(l + u) uv(u + v + uv)—u(S x uvS 1 +uS l +S~ 
2 \v(SuvS + vS + Su)u-uv(u + v + uv) vSuv(\ + v) - uv(l + v)S~ l v 

fi = In m or ji = In v 
• Recursion operator: 

_j / uv-\S~ l + u + v + uv u(l + u_i)5 _1 + m(1 + u) 

K = <H 2 <H l =\ v[l + Vi)S+ u^v A u\vS + ii\ + v + U\V\ 



1 v)v\ 



U(V - V-\ +UV - U-iV-i) 
V(U\ -11 + U\V\ - uv) 



Non-trivial symmetry: 

( UV(1 + U)(U + U\+ U\V\) + UV 2 (1 + U) 2 - UV^Jl + U-i) 2 - U 2 V-i } 
Ut \ _ -U-\UV + U + V-2 + U-i + U-2V-2); 

Vt J ~~ U\VV\(\ + 2ll\ + + Ul) + U 2 V + U\VV 2 (\ + U\) + V 2 U\(\ + V\) 

-UV(1 + V)(V + V-\ + U-1V-1) - u 2 v(l + v) 2 



n 



Lax representation 



M 
U-- 



A 2 + 2ui +2v + l -A(2v + 1) - A~ 1 (2u 1 + 1) 
A{2u l + 1) - A~ 1 (2v + 1) A" 2 + 2«i + 2v + 1 

A 8 A + UV + \ + : 
1(20+1) 



A^ + A 



A £E±i) +A -ieH±I) 



A 2 -A- 



+ UV + j + j 



In [65J, the Lax representation is given in the following form: 

U t = U C - A U; W t = WB-CW; V t = VB-AV, 

where the difference operators U, V, W, A, B and C are given be 

U = u + AS~ 1 ; W = l + A~ 1 vS; V = 1- A~ l uvS 

A = u + U-1V-1 + V-\ + AS~ l ; B = u + uv + V-\ + AS~ X ; C = u + uv + v + A<S _1 
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There exists another weakly nonlocal recursion operator The recursion operator % has a weakly nonlocal 
inverse: 



( 



uv o , i»'("+l)+((l+t'-i) 2 +")(»i+l) 



»(»+!) 
' (ll+V-i +1) 2 ' 



-1 



5- 



M+Z7-1+1 Ml +0+1 
«+!)_!+! Ul+V+1 



u(u+U_i+l)(Ui+D+l) 

(S-l)-^ 



2 1 



2 



»(»!+!) 

(i+m)(i+p) 



(m+H+i) 2 



In fact, it commutes with H and satisfies ??' + id) = id. 

It is related to the Relativistic Toda equation l|83t , written in variables u and V, by the Miura transformation 
u = -uv and V = —(u + V-\ + 1) [66|. This transformation is similar to j86l l, which explains the name of 
this equation. 



4.10 The Merola-Ragnisco-Tu Lattice 

• Equation IT51I671 : 



U t = Ui - ifv 



• Hamiltonian structure JT5): 



1 
-1 



-V-\ + V u 



f = mv- 



1 ? 

irv A 



• Recursion operator fl5l : 

n = 

• Non- trivial symmetry 11151 : 



S-2uv 



n 



Master symmetry lfl5l : 



Lax representation: 



u t 

Vi 



(n + l)u 
-nv 



M = 



+ 2 ^ \(S-l)-*(v u) 



ii2 - u\v\ - u 2 v-i - 2umi\ + u 3 v 2 
>-2 + v\u.-i + V 2 U\ + 2lWV-i - U 2 l 



nu t + 2ui - 2u 2 v - 2u(S - l)~hw 
nvt + v-i + uv 2 + 2v(S - Yf l uv 



v 

-2A - uv 



U = 



-A 

-u 



-v-i 
A 



The recursion operator H has a weakly nonlocal inverse: 



i 



(lLjD+1) 2 



+2 



(uui+1) 2 



U-\V+\ 

__Ei_ 

UVi+1 



(ud + 1) 2 



(H_ 1D +1) 2 



The symmetry (u, -u) tr is a seed for both H and 7? 

Under the invertible transformation 1 i— » —t, u i— > — M and ph d 1( this lattice transforms into 

M f = -Mi - U 2 »i 
= C_i + 2? 2 M_i ' 

which is related to the nonlinear Schrodinger system u t = u xx + 2u 2 v, -v t = v xx + 2v 2 u appeared in l68l . 
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4.11 The Kaup Lattice 

• Equation 



• Hamiltonian structure 



Ut — (u + v)(u\ - u) 
v t = (u + v)(v - V-l) 



u + v . 
1 -\u + v) 



• Recursion operator: 

I (U + V)S + U 1 -U \ I U t L_ ,4/ j j \ 

+ | \ ^S(S-VT X ( v-i-v ui-u) 

• Non-trivial symmetry: 

<fcl U > \ _ / ( M + v )( uu l + uc, -l + u l v l ~ u l u 2 - «2^1 ~ U\V-\) 
\ V t J ~ I (U + v)(U-iV-2 + V-2V-1 - U\V - U-1V-1 + U\V-\ - V-\V) 

• Lax representation l68l : 



1 v-A ' 1 v-i 



There exists another weakly nonlocal recursion operator 

( U+V O-l (»-»-! ) 



ft' 



(u-i+vf (u-i+v) 2 

(Pi -p) il+P O . »-»-i-Pi+P 

k (h+pi) 2 (m+^j) 20 (a_i+o)(«+t>iV 



/ 1 \fs_ n-ifj L. _J L.^ 

' \u+Pi U+P K-l+B H+P^ 



»-»-! \ 

u-i+o /<?_n-l/J L_ _2 

IV-' / \M+Pi U+P U_!+P U+V ) ■ 

U+l'l I 



The symmetry (1, -1) is a seed for both K and 7?'. In fact, operator JZ' is the inverse operator of 
4.12 The Ablowitz-Ladik Lattice 

• Equation 

' ut = a-uv)(au 1 -(lu- 1 ) :=aR 
Vt = (1 - uv)(pv\ - av-i) r 



• Hamiltonian structure llT5l : 



Recursion operator IT51I7D1I7T1 : 

_ / (1 - MP) 5 - Mi t> - WU-i —UU\ 1 / -K Ln <v-1/ \ 
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• Non- trivial symmetry Ifl5l : 

(1 - uv)u\ 
-(1 - uv)v-\ 



7? 



(1 - UV) ((1 - U\V\)U2 - VU 2 - UUiV-ij 
(1 - UV) (-(1 - U-iV-\)V-2 + UV 2 ^ + U\V-\v} 



Master symmetry 1IT51 : 
<R 



nu \_l (n + 1)(1 - uv)u\ - u 2 V-\ - u(S - 1) x uv-\ 
-nv J I (1 - «)(1 - uv)v-\ + uvv-i + c(»S - l)~ l uv-\ 



• Lax representation: 



. A u \ 1T ( A 2 - uv-i Am \ n I A 1 m_i 



The coefficients for a and j8, namely Xi and K-\, are commuting symmetries for the equation. The 
inverse of this recursion operator "R is of weakly nonlocal form: 

\^ -VV\ (1 - MD) <> - WZ?j - U-iV ) \ -v ) y ' V > 

/ (l-uv)u-r \,o_ ir iC _«_ \ 

\ -(1 - hz;)i>i j ^ J V l-w l-«w J 

Both 7? and share a common the seed cr = | j' Starting f° rm it, we can generate the commuting 

symmetries 7?~'(cr) and 7?'(cr) for z e N. For more Lie algebra structure among symmetries and master 
symmetries, we refer the reader to [15]. 

4.13 The Bruschi-Ragnisco Lattice 

• Equation 1721 1551 l40l: 



ll t - U\V — UV-\ 

v t = v(v - V-l) 



Hamiltonian structure l72l[T5l : 



\v{S-l) )' ^ = U ' V 

i vSu - uS~ l v v(S -l)v \ . 
^ = \ v(l-^)v ' f 2 = u 



Recursion operator IliOlUBI : 

K - " 2 -\ vS- 1 + v(v - v-i)(S - I)" 1 1 



• Non- trivial symmetry lfT5l : 

n 



U t \ _ I VViUj - V-2V-1U 
V t J ~ I V(VV-1 - V-1V-2) 
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Master symmetry ITBI : 



n 



nu t + 1u\v + U-\u 
nvt + vv-i 



The recursion operator H has a weakly nonlocal inverse: 



C-l 1 C-l U 11 i UVO -I \— 1 1 

(D-\ _ qj siJ-l _ | ° p ° S 2 " ^ + ^~ p A»> ^ p 







_ c-l _u u_ \ / IM. _ u \ 



However, recursion operators H and H 1 have different seeds similar to the Relativistic Toda system in 



section l47l The starting point for operator H is the equation itself while the seed for R l \so = y "v l \ y 

Operator fi~ l acting on the equation itself and "K on a do not generate new symmetries. We refer the 
reader to [15] for more Lie algebra structure among symmetries and master symmetries. 
In fact, the Bruschi-Ragnisco lattice is trivially solvable, so are its symmetry flows. The equation for the 
second component t>t is independent of the first component u. Besides, the scalar lattice v t = v(v - v _i) 
and Vt = ^- — 1 can be linearised into zv t = W-\ and zv t = W\, respectively by the transformations v = — ^ 
and v = respectively. Once it is solved, the equation for u is then linear. 



4.14 The Kaup-Newell lattice 



• Equation Il73l : 



= a(— ^ —) + b(— 



:= flKi + 



• Hamiltonian structure l73l : 

/ S — 1 \ 
H - I q ' / = -aln(l - wc) + frln(l + wci) 

• Recursion operator: 



<R = <HI = 



i 



(l-UlHl) : 



1 



2»iP 



(1-UP) 2 

-ius-i)-^ rfs ^ ), 

where the symplectic operator J is defined by 



(l-iro) 2 (1-HiPi)(1-up) (l-idm) 2 ^ T (1-iro) 2 (l-Kp)(l-iiiPi) 
-5-1 _ -Jt_ 1 C-l + JL=2««L 



(1-i^P-iF 



(1-lfp) 2 



I = 

Non-trivial symmetry : 





l 

l—uv 



1 

l—uv 





1-uo 



^! (5 + i)(5-ir 1 ( ^ ^ ) 



WO 



1 /„ _ »2 "l P \ _ 1 /,, _ »1 » 2 P-1 \ ^ 

(1-HlPl) 2 1-H 2 P 2 (1-UU) 2, > I-HjCj l-H_lP_l ' 

J /„ »lP 2 P-l \ 1 ""-1 P-2 \ 

k (l-l(p) 2 » 1— «iPi l-H_iP_i' (1— H-jP-i) 2 * ~1 1-UP 1— U-zV-2' , 



• Lax representation (73): 

A + (1 - A)mp u 



M = 



U = a 



(1 - A)v 



A-l 
2 



1— MP 



(l-A)p-i A-l 

1— H_i P_i 



(A - 1)(1 - H_ lP ) 



2m_ 



^ / 2A(1 + m_io)1 2(1 -A)p (1 - A)(l - tt-io) 
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The recursion operator % has a seed a — | v \ anC ^ = ^ ^ a ^ n § integration constant being zero. 

Similar to the Ablowitz-Ladik Lattice in section 14.121 the coefficients for a and b, namely, K\ and K-\, 
are commuting symmetries for the equation. Indeed, there exists another weakly nonlocal recursion 
operator 



<R' = <HT = 



' 1 £-1 _ l+2»Pi u 2 £ + "-i 2»»_ 



(l+H-lO) 2 (l + «Ol) 2 (l+!(fl) 2 (1+H-lt?) 2 (1+H_lt»)(l+H0l) 

p 2 ^-1 _ "i 1 1 2m_iPi 



V (1+U-lt)) 2 (l+«»l) 2 (1 + HOi) 2 (1+U-lV) 2 (l+M_lO)(l+HBl) / 

-2K_!(5- 1) _1 ( ), 
where the symplectic operator I' is 



v^S s - u -^)uk-v 



1+u 



Again a is a seed for 'R' and 7?' (a) = K_i taking 

(5 _ irl /_«5_ + «^_\ 

V 1 + UV\ 1 + U-iVI 1 + U-\V 

By direct calculation, we have I"R = IH' - I' — I. Thus these two recursion operators satisfy the 
relation K"R = KR! =K -% that is, {K - id)(9? - id) = (9? - id)(9?' - id) = id. 

4.15 The Chen-Lee-Liu lattice 

• Equation 1751: 

u t =a(l + uv)(u 1 -u) + b(l + u- 1 v)- 1 (u-u- 1 ) 

Z7f = + - V-\) + b(l + UV\) L \V\ - V) 

• Hamiltonian structure: 

1 + uv 



(U I 1+UV \ , . X71 1 + 

hi = I > „ , / = a(uv-i - uv) + b In - — 

^ -(1 + ut;) /' 7 v ' l + 

• Recursion operator: 

<n <u I (1 + Mf)«S - 2«C + U\V + MZ?_i UUj - « 2 \ 

\^ v z - VV-i (1 + uv)S 1 j 

where Hamiltonian operator ^2 is given by 

/ (1 + mu) (5(1 + uv) - UV + U\V) \ 

^ 2 = \ (uv - iiiv - (1 + uv)S~ 1 ) (1 + uv) j 

• Non-trivial symmetry : 

^ / (1 + UV){U\UiV\ + Uz + M 2 l> + UU\V-\ + U 2 V - U\ - « 2 U_i - 2llU\V - U 2 V\) 

1 (1 + UV)(ll\V 2 + 1uV-\V + U-\V 2 _ X + C_i - C_2 - Wf 2 - iro 2 x - U\V-\V - U-\V-2V-\) 
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• Lax representation |73| : 

A + uv u 

M = 



(1-A)v 1 ' 



U = a 



A - 1 + uv-i u 



(1 - A)t>_i A(l + m-io) \ (1 - A)o 1 - A 



The coefficients for a and £>, namely, K\ and are commuting symmetries for the equation. The above 



recursion operator H has a seed a = 



and K{o)=K l . 



There exists another weakly nonlocal recursion operator 



<R> = <H'<H- X = 



l+uv c-1 



K-l (»-»-! ) 

" (1+U-ip) 2 

Pl(Pl-P) 1+Mg C , Pl»-2»-iPi+»-iP 

'(1+UPj) 2 (l+iiPj) 20 (1+H-xOXl+HPl)^ 

2lL 



ix ( o -1 \-l / 2pi p_ 2h_i » \ 

-l _M W«_ ir 1 ^- 3 v — M 

I j ^ A / ^l+UPi 1+UP I+U-1P l+uvj' 



where the Hamiltonian operator fi^ is 



<H' 2 



n l+»p / C-l 1+KP , \ 

1+H-lP V° 1+UPi 1 + UPi J 



_ ( 1+»p C . Pi(»-»-i) \ 1+up 

V \l+MPi 1 + HPj ) 1+U-lV 

+ 1)(5 - l)" 1 ** - X_ 1( S(5 - l)" 1 ( -« o ) - 1 ~" J (S - l)" 1 ** 
Again a is a seed for H' . In fact, operator H' is the inverse operator of H. 

4.16 The Ablowitz-Ramani-Segur (Gerdjikov-Ivanov) lattice 

• Equation 1Z3: 

Wf = (a«i - + mz>)(1 - wci) 



v t = (bv\ - av-\)(l + uv)(l - u-\v) 



:= aKi + bK. 



Symplectic operator: 



I = 





_j p-i l 

1+MP \—UV\ 



i g L_ 

1-i/Pi 1+WP 







and we have 

I(aK\ + bK-i) = 5( UfV ) (a(uv-i - uv - ux>U\V\) + b{u-\V\ - uv\ + U-\uvvi)) 
• Hamiltonian structure: 

\ -(i - + wujo ' + uv-i)(i + uv) 



(1 + MZ>)(»S(1 + MZ>) + UV-i)(l - U-\V) 

(1 - u-\v)((l + uv)S~^ + uv-i){l + uv) 



-KrS(S -!)-'( u -v)-l \ J^-l)- 1 ^ 
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Recursion operator: 

' (1 + UV)(1 - UVi)S + U\V - U\Vi -UUi(l + UV)S- M 2 (l + U-iV-i) s 

+uv-i - uv(l + U-iV-i + 2u\V\) + 1 1 T„" Pi p mi(m - u-i - 2uu-\v) 



n = 



+ 



(1 - U-\V)V-\V - (1 + UV)V-\VS 1 (1 + UV)(1 - U-iV)S 1 + UVU-iV-i J 

~ uv \) _J> £i_ _« »-i \ 

— U \V) J V 1+uv 1-a-jD J 

- C_i + U-\V-\V + U\VV\ U -U\+ UU-\V-\ + UU\V\ ) 



«i(l + uv)(l - UV\) 
-V-i(l + uv)(l 



= 911 + 



1 



1 

• Non- trivial symmetry 

((1 + UZ?)(1 - UVi)((l + Mi?7i)(l - U\V%)Ui - u\v\{l + iro) + UU\V-\{1 - U-\V) - U\V(ll - Ml)V 
(l + M?7)(M-lZ>-l) ^(l + M_ii;_ 1 )(l-U_2!7-i)^-2-M-l^ 1 (l + Ml') + MlU-l^(l-iroi) + MU-l(u_i-^)^ 

• Lax representation 1731 : 
M 

\ I 

jj _ a ( & + M?3-i(l - 



The equation given in (73l is 



A + uv u 
(1 - A)(l - uv\)v \ — ux>\ )' 



(1 - A)(l - u-\v)v-\ -uv 

b I ii-\v - A u-\ 
A \ (1 - A)(l - ii-\v)v (1 - A)(l - u-\v) - \u-\Vi{l + uv) 



Hi 



(b - a)u 
{a - b)v 



+ aK\ + bK-i. 



Since the vector a = commutes with both K\ and K-\, we removed this term in our consideration. 

\ v ) 

There exists another weakly nonlocal recursion operator 



H' = 



(1 + UV)(1 - UV\)S 1 + UVU-\V\ 



UU-i(l + UV)S + UU-2{\ + U-iV-i) 



1-UVi 



(1 - U-\V)V\V + (1 + uv)v\vS 1 
+ UV){1 - UV\) 

V\{\ + uv)(l 



'-^u-i{u-u-i - 2uu-iv) 
(1 + uv)(l - u-iv)S + uvi - 2u-\uvv\ 

-U-\V\ + U-\V - U-2V{1 + U-\V-i) 

— uv\) —2 — a "-1 \ 

U \V) J V l +uw X—UV\ 1+uv 1—U-\V J 



l -"p j (*5 — 1) 1 ( ^2(1 + U\V\) + V\(U-\V - 1) W_2(l + W-lC-l) + U-i(uVi - 1) ) 



where the Hamiltonian operator W is 

<H' 



-(1 + uv)(l - u-\v) 

(1 + uv)(l - u-\v) 



-K^SiS u -v)-l ^ 
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Operator 9?' is the inverse operator of %. The vector a is them seed for both of them and K'iK-i) is 



((1 - UV\)(\ + UV) [{1 + U-\V-\)(1 - U-\V)U-2 - UU-\V2{1 + U\V\) - U^Vi(l + Uv) + U-\(U-\V + UVijj 
(1 - U-\V)(1 + Uv) (-(1 + U\V\){1 - UV\)V2 + U-2VV\(1 + U-\V-\) + U-\V\{V\VU + V\ — v) — uvfj 



4.17 The Heisenberg f erromagnet lattice 



• Equation IITil : 



Ut = (m - - Mi)(mi - 1 

C t = (M - - V)(U - 



• Hamiltonian structure: 

1ri 



(M-57) 2 . 

-(u- v f o I' / = ln(u - o) - ln(u - 



Symplectic operator: 



I 



/,. _ 7 ,\-2 C _ 
^ W V) O (u _ D) 2 (u _ ; ,_ l)( „ 1 _ I , ) 



)(»l-f) 







»-»! V'-' + -"-A I («-t;)(u-P-l) («-»)(«! -») ) 

(H-B)(«l-C) 



Recursion operator: 



( (u-v) 1 



<R = <HI = 



S- 



2(»-»i)(p-P-i) 



(»-"i) 2 

(Hl-f) 2 

(k-p) 2 q-1 



(«i-h) 2 (u-h_i)(«i-d) 

(p-p-l) 2 

j (vS - 1) 1 ( (u-^^-Vj) („-")("} ) 



Non-trivial symmetry: 



9? 



/ (u-v) I (U-V)(u x -Vi)(ui-U2) (u-Uif(V-i-V)\ \ 



nr(v-i-v) \ 



(Ui-I)) 2 V ("2-»l) 

(»~P) / (»-P)(»-l-P-l)(P-2-P-l) + (P-P-l) 2 (»~»l) \ 
, (M-D_l) 2 V (U-l-P-2) («l-f) / y 



• Lax representation: 



M = 



A - 2u(u - v)~ l -2(u - v)~ l 
2uv(u - v)~ l A + 2v(u - v)' 1 1 

i) 



-,-1 u + v-i 2 
v ; \^ -2uv-x -(u + v-i) 



The recursion operator H has a weakly nonlocal inverse: 



K' 1 = <HT = 



( (u-vf Q-l 

(»-i-«) 2 ° 



(«-»-i) 2 

(u^-v) 2 
(u-v) 2 q 2(w— H_i)(P— V\) 

(ll-Vi^U-x-V) J 



\ (u-ih) 2 (u-vi) 2 
f (u—v)(u_\ —ll) \ 

^° J-- 1 ^ (U-V)(u-Vl) (U-V)(U-!-V) )' 



U-i-V 

(u-v)(v-vi) 



\ u-vt 
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where the symplectic operator I' is given by 



r = 



o -(u - v!)- 2 s + . ( ".;,"- i)(; ;: Pi) . ^ 

(u_i - v)- 2 S^ - , ("-"-'X''-" 1 ) 

v 1 ' (u-vy-(u— V\)(U-X— v) 



(87) 



H«-> + i;(,o i; ^ (u-v)(u-v{) tu-v)(u- 1 -v) ) 

(u-v)(u-i-v) I 

The Heisenberg ferromagnet lattice is a special case of the following Landau-Lifshitz (Sklyanin) chain 

where 

fc( M , p) = I(JCj(l - uvf - K 2 (l + «u) 2 + K 3 (u + vf), K!K 2 K 3 #0, K„eR 

by specifying a = l,b = and /i(w, p) = i(u — v) 2 . Notice that the coefficients of a and b in the Landau- 
Lifshitz (Sklyanin) chain j87t commute. 
The Lax representation for equation d87t is given in l74l by 



!J = i Y j aN^(A)S k (u,v- 1 )a k + bN"(A)S fc (i; / M-i)cT fc ; M 
Here 0^ are the Pauli matrices 

0-1 = ( 1 )' ° 2 
vector function S(p, q) is defined as following 



-i 

1 /' 



V(SKS) 



(7 3 



Y_ l M k {\)S k {u,v)a k 



I 

V k=i 



1 
-1 /' 



(S 1 (p,q),S 2 (p,q),S 3 (p,q)) = 



1-pq i + ipq p + t 



p-q p-q p-q 



satisfying S 2 + S 2 + S 2 = 1, 

3 

(SKS> = 2^K ; S 2 (w,z>) 

and M;(A), Nf(A) can be expressed in terms of Jacobi elliptic functions of a spectral parameter A: 
Mj (A) = Jl-KiK-hni^x), M 2 (A) = JT^^cn(A, k), M 3 (A) = Jl - K 3 K 1 _1 dn(A / k) 



N 1 ± (A) = ^S Ml (A± H ), NJ(A) = Sm 2 (A ±/J ), N*(A) = 



M 3 (A + /.<), 



where ?c and /.i are defined by equations 



K 3 (Ki - K 2 ) Ki 
K= V^lj' ™^=K 3 

Instead of explicit unif ormisation of the elliptic curve, one can use identities 



M 2 - 1 M; - 1 



N,± = 2(M 2 -1) (M ' * M/M °' ? ' >' k G {1, 2 ' 3} ' l * ] * k + 1 
One local Hamiltonian operator for equation d87t l74l is 

1 



t H = h(u,v) 



-1 ' 



v) h(u, v) 

f = a In -r- 1 7T + b In 



(«i - p) 2 



(M-l-f) 2 ' 
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Knowing the Lax representation, we can, in principle, compute its recursion operators (as we did in 
section I3T3 with multipliers 

M A ) = (Nf~(A)) 2 and v ± (A) = (A)Nf (A)N*(A). 

However, the calculations involved are rather big and we haven't found a neat way to present the 
operators. 

4.18 The Belov-Chaltikian lattice 

• Equation [76): 



U t - U(V2 - V-\) 

V t = ll-i - U + v(t)\ - V-i) 



• Hamiltonian structure Il76l : 



= lu{S-S- l ){S + l+S- l )u u{S - 1)(S + 1 + S-> 
1 " I 0(1 - 5~ a )(5 + 1 + S'^u viS-S-^v + S^u-uSj' 



fi = v 



<H 2 







i(l + S + S 2 )(uS - S~ 2 u) 



(uS 2 - S^u) (1 + S' 1 + vS" 2 )u o(l + S)(uS - S~ 2 u) + (uS 2 -S-hi)(l + S'^v 
u(l +S + S 2 )(S-\>-vS)(l + S' 1 + S'^u u(l +S + S 2 )(S-\!-vS)(l + S~^v\ 



v(l + S)(S-h-vS)(l + S- 1 + 5-> 



v(l + S)(S-h>-vS)(l + S-^v 



h = "3 ln " 



Non-trivial symmetry : 

<H 5f =( UV ~ l (° + U_1 + ^ ~ "° 2 ^ +V 2 + V 3> + U ("1 + "2 - M-l - U- 2 ) 
2 * 1 ~~ 1 (tt - t7l7j) (Z7 + 17} + 172 ) + (t7t7_i - (l7 + t7_i + Z7_ 2 ) - 17 («_ 2 - Ml) 

Master symmetry 11771 : 

nut + uv\ + 4ui7_! + wz7 
nvt + u — vv\ - 4 u_i + 4 i7i7_j + v 2 



• Lax representation (78ll : 



' A Ai7 A«_i 



M = 



( v - A 



U = 



-Av -Au-\ 
v-x 

-1 Z7_2 



The Belov-Chaltikian lattice is the Boussinesq lattice related to the lattice W3-algebra (78l . 

In recent paper [79 1, the authors wrote down the Boussinesq lattice related to the lattice W,„-algebra for 

the dependent variables u , u 2 , ■ ■ ■ , u m and independent variables n and f as follows: 



u) = -u\u 2 m - u 2 _ x ) 
' -u< 
u m (u 

-1 

The vector t = (t 1 , • • • , t") t defined by 



u\ = u m - u l +l - u l (u 2 _ x - u\) r i = 2, 3 • • • ,m - 1 



2 

m-1 



t l =nu\-u l ((m + l)ul + Zl u 2 ); 
t' = nu\ - u' (iuf^ + Y!i=o u 2 ) + (i + l)i 



i = 2, 3 • • • , m - 1 



x m = nu n t - u m {mu 2 m _ x + u ]) + ( m + i)" 1 
is a master symmetry. Its Hamiltonian structures are also studied in the same paper. 



35 



4.19 The Blaszak-Marciniak lattice 

• Equation ED: 



U t = W\ — W-\ 

V t = U-ilV-i - uw 

iff = w(v - V\) 



• Hamiltonian structure 

( s-s- 1 










(S- 1 - l)w 




<H 2 



-w{S - 1) 

'Sv-vS^-uiS + ly^l-Sfu SioS-S-hv 
= wS - S~ x wS~ x 

w(S + l)- l (l-S)u 

f 2 =v 



fl = uw + -v 2 



S l uw—uwS 



u(S + l)- 1 ^ - S)w 
viS' 1 - l)w 



-w(S - l)v wiS^-Sjw-wiS + 1)-\1-S)w 



• Recursion operator: 
where 

• Non-trivial symmetry : 

'Hibfi = 'Hi 

• Master symmetry: 



I 1 










o -{s-\y l i 



( w ^ 

v 

y u 



lfl(27i + l>2) ~ W-\(V + V-l) 
U-\W-\{V + V-i) - Uw(v + V\) - W-1W-2 + WVJ\ 
IV(V 2 - V 1 ) + ZV(W-iU-i - W\U\) 



u/2 v 
v 

I 3w/2 J 



• Lax representation 1801 : 

L = S 2 + u\S — v\ + wS~ l , A = S 2 + uiS - v\ 

We do not explicitly write out its recursion operator, which is no longer weakly nonlocal although both 
operators 'Hi and < 7Y~ 1 are weakly nonlocal. The statement that such recursion operator generates local 
symmetries can be proved in the same way as in l30l , for weakly nonlocal differential recursion operators. 
We can compute the next Hamiltonian is uvw + uv\W + y - ww\ . Its master symmetry is highly nonlocal, 
which is explicitly given in (77ll . 
Another three-component lattice was given in Il28l as 

p t = q 1 -q 

q t = q(p-i -p) + r-r- 1 . 
r t = r{p-i - p-i) 



It has the Lax representation 



M = 



( 1 
q u + A 1 



U = 



' -p-i 
<? 

r-i 



1 

A 




, r 

This lattice is related to the Blaszak-Marciniak lattice by the Miura transformation 

p = v\, q = -uw, r = -ww\. 
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5 Conclusion 



In this paper, we review two close concepts directly related to the Lax representations: Darboux transfor- 
mations and Recursion operators for integrable systems. We use the well-known nonlinear Schrodinger 
equation, whose Lax representation is polynomial in the spectral parameter, and a deformation of the 
derivative nonlinear Schrodinger equation, whose Lax representation is invariant under the dihedral 
reduction group D2, as two typical examples. We then present a list of integrable differential-difference 
equations containing equations themselves, Hamiltonian structures, recursion operators, a nontrivial 
generalised symmetries and Lax representations. For most equations, we also add notes on their rela- 
tions with other known equations and the weakly nonlocal inverse recursion operators if existing. 
The theory of integrable partial difference (or discrete) equations is a relatively recent, but very active 
area of research. There are famous the ABS equations for affine-linear quadrilateral equations due to 
Adler, Bobenko and Suris I5TII52"! , based on the condition of 3-D consistency. Then Levi and Yamilov 
proposed to use the existence of a generalized symmetry as a criteria for integrability [83] to classify 
integrable partial difference equations. They extended the ABS list. For all integrable partial difference 
equations, their symmetry flows can be viewed as integrable differential-difference equations. The 
concept of recursion operator has been adjusted to the difference equations to show that it generates 
an infinite sequence of symmetries and canonical conservation laws for a partial difference equation 
l34l [Til l52l . It can be proven that the difference equation shares the same recursion operator for its 
symmetry flows. Thus the study the symmetry structure for integrable difference equations is in fact the 
same as for integrable differential-difference equations. This is one of the motivations to produce the list 
presented in this paper. Besides, this list can serve as a benchmark for developing computer software 
packages for the symbolic computation of recursion operators, Lax representations, symmetries and 
conservation laws for nonlinear differential-difference equations f7T\ . 

We have not presented here a rigorous proof that the operators we computed from Darboux transfor- 
mations are really Nijenhuis recursion operators for the equations obtained from the corresponding Lax 
representations. When the Lax representations are polynomial in the spectral parameter and under 
certain technical conditions a sketch of the proof was given in 1 28] . We think that the main statement of 
EHl can be greatly generalised and simplified and there is a neat rigorous algebraic proof that the op- 
erators obtained from Darboux transformations invariant under reduction groups are indeed Nijenhuis 
recursion operators. 
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